Anomalous low temperature ambipolar diffusion and Einstein relation 
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Regular Einstein relation, connecting the coefficient of ambipolar diffusion and the Dember field 
with mobilities, is generalized for the case of interacting electron-hole plasma. The calculations are 
presented for a non-degenerate plasma injected by light in semiconductors of silicon and germanium 
type. The Debye-Huckel correlation and the Wigner-Seitz exchange terms are considered. The 
corrections to the mobilities of carriers due to difference between average and acting electric fields 
within the electron-hole plasma is taken into account. The deviation of the generalized relation 
from the regular Einstein relation is pronounced at low temperatures and can explain anomaly of 
the coefficient of ambipolar diffusion, recently discovered experimentally. 

PACS numbers: 71.27.-|-a,71.35.Ee, 78.20.Jq, 78.56.-a 



Introduction 

The Einstein relation for electrons Jj nUdji/dn = eD 
connects mobility U with diffusion coefRcient D by the 
thermodynamic function dfi/dn and electron charge e. 
Here fi is chemical potential, n is electron density. For 
the case of the Boltzmann gas of non-interacting elec- 
trons, where dfi/dn = kT /n, one gets a regular form of 
the Einstein relation, presented by Einstein [2l] and von 
Smoluchowsky Q for the Brownian particles. Here kT 
is the temperature in energy units. However, the Ein- 
stein relation can be used in a general case of interacting 
particles[l]. 

The focus of this paper is the ambipolar diffusivity 
under condition that interaction between electrons and 
holes is substantial. I have been initiated by the paper 
of Hui Zhao^ who has measured the coefficient of the 
ambipolar diffusivity (CAD) by optical method in silicon- 
on-insulator (SOI) structure with 750 nm silicon layer. 
The doping density is lO^^cm"'^ while the density n of 
electron-hole pairs excited by light is between (0.5 — 3) x 
lO^^cm"^. The temperature is in the range ^OK — AQOK . 

The regular Einstein relation for the CAD in non- 
degenerate and non-interacting electron-hole plasma has 
a form 



2kT UpUn 

e U„ + U„ 



(1) 



where C/„, Up are mobilities of electrons and holes respec- 
tively. The CAD observed in Ref.[il approximately fol- 
lows Eq.lID at 400i^ > T > 300K. At lower tempera- 
tures the CAD goes down and deviates from the Eq.([T]) 
approximately 7 times at T = 90K. The mobilities were 
taken from experiments with a pure bulk silicon. Author 
explained this anomalous behavior as a result of uncon- 
trolled defects that are absent in a bulk silicon. 

I propose here alternative and more universal intrin- 
sic explanation of the low temperature anomaly. This 
explanation is based upon a novel form of Einstein's fun- 
damental relation for the CAD, that takes into account 
interaction between excited carriers. 



Under conditions of the experiment the electron-hole 
plasma can be considered as a non-degenerate. On the 
other hand, the classical correlation energy per particle 
^ g2^i/3^^ jg ^jjg order of 200K. Here k is dielectric 
constant. Thus, in the region of the observed anomaly 
the interaction energy becomes of the order of tempera- 
ture. 

The paper is organized as follows. First the novel form 
of the fundamental Einstein relation connecting the CAD 
and the Dember field with mobilities of electrons and 
holes is derived. The relation contains the derivatives of 
the Helmholtz energy density (HED) of the interacting 
plasma with respect to particle densities. The HED is 
calculated taken into account correlation and exchange 
between particles. Then the corrections to the mobilities 
due to deviation of acting electric field from the applied 
field are considered. Finally the theoretical results are 
compared with the experimental data. 



Einstein relation for ambipolar diffusivity. 
Thermodynamics approach 

Silicon is an example of semiconductor with a long re- 
combination time of interband excitation. Assume that 
this system is in the thermodynamic equilibrium with re- 
spect to all relevant parameters except the total numbers 
of electrons and holes. To derive the Einstein relation I 
use here the same method as in Ref([4]). The difference 
is that excited electrons and holes have two independent 
electrochemical potentials $" and $p respectively. The 
Helmholtz energy has a form 

P = J f{n,p)d^r + J e(p(r) - n(r))V'dV 



$P / p{r)d^r - $" / n{r)d 



(2) 



Here n and p are electron and hole densities. Function 
f{n,p) is the HED of the almost neutral and microscopi- 
cally homogenous electron- hole plasma, the function V'(r) 
is a potential of a static electric field. 
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Using conditions 5F/6p = and SF/Sn = one gets 



and 



op 



on 



(3) 



(4) 



It follows from the general principles of statistical physics 
that that in the state of equilibrium both <i>P, $" should 
be constant along the system. Exploring Einstein's idea 
that electric field is equivalent to a certain density gradi- 
ent one can write the fluxes of holes and electrons as 



-^V<i>^q„ = -^V<I>" 



(5) 



where dp and an are conductivities of holes and electrons 
respectively. Then the fluxes are 



and 



E 



dp^ 



9V 



dpdn 



dpdn 



Vp 



(6) 



(7) 



Note that separation of the conductivities of electrons 
and holes are possible only if their mutual scattering is 
small. We assume here that the mobilities of the carriers 
are controlled by the lattice scattering. But even in this 
case the above expressions predict a drag effect due to 
the interaction terms in the HED. The flux of holes is 
proportional to gradient of electron density and to gra- 
dient of hole density. The same is true for the flux of 
electrons. For example, if electric field is zero and gra- 
dient of electron density is zero, there is an electron flux 
proportional to a gradient of hole density. 

The ambipolar diffusion is measured under condition 
that electrical circuit is open. Then the total electric 
current j is zero and Qp = Qn- Taking into account the 
continuity equation 



ed{p — n) jdt -I- divj = 0. 



(8) 



one finds that since electron and holes are excited by 
light in equal amounts the system is neutral everywhere. 
Then n(r) = p(v) and Vn — Vp. A small separation of 
charges at the boundaries of the sample appears due to a 
difference of diffusion coefficients of electrons and holes. 
This difference is compensated by an electric field called 
the Dember field. Nevertheless, the electron-hole plasma 
in the bulk of the sample is neutral. 

The expressions for flaxes can be written in a form 



qp = ^E - DpVp 



(9) 



and 



— E - A,Vn. 



(10) 



It follows from Eas. (|6l7l9ll0p that diffusion coefficients 
of electrons and holes Dn and Dp are connected with 
corresponding mobilities t/„ , Up by the relations 



and 



" \ dn^ dndp J e 



^ \ dndp J e 



(11) 



(12) 



Eqs. (lllll2p give a novel form of the Einstein relation that 
is valid in the case of interacting plasma. 

The mixed partial derivatives in these equations de- 
scribe the drag effect. 

Using Eqs. ()9|10p and condition qp = q„ one gets ex- 
pression for the Dember field in terms of D„, Dp 



E. ^ <£l^Vn. 



(13) 



This regular form becomes more complicated if Dn and 
Dp are expressed through mobilities C/„ and Up using 
Eqs.iniini). Then 



E, 



[e{Un + Up)Y 



1 W 
dn^ 



Un 



u„ 



dp^ " 



dndp 



[Un - C/p) V 



(14) 



Due to the interaction of carriers the Dember field is 
not necessarily proportional to the difference of mobilities 

Un - Up. 



Substituting Eq.dTS]) into Eqs. ((9TT0| one gets 

Qn = Qp = -DaVn, (15) 

where the CAD 

DnUp + DpUn 



Da^ 



Up + Un 



(16) 



Using Eqs. (|llll2[) we get the generalized Einstein relation 
for the CAD 



2kT_UpU^ 

Da = ' Q{n,T), 

e Up + Un 



where 



Q(n,r) = 



n fd^f , d\f ^ ^d^f] 



2kT \dn'^ dp^ dndp J 



(17) 



(18) 



is a ratio of the coefficients of ambipolar diffusivity cal- 
culated with and without interaction (cp Eq. ()17|) with 
Eq.©).^ 

It is important to put n = p after calculation of the 
second partial derivatives in Eqs. fTTl 112114118]) . 
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HED of semiconductors with band structure of Si 
and Ge 

Analytical calculations of the HED of interacting carri- 
ers are possible in the framework of perturbation theory 
only. One should keep in mind, however, that the region 
of applicability of these calculations does not cover all 
temperature range of the experiment. 

The HED can be written in a form / = fid + fi, where 
the first term describes the ideal gas, while the second 
one takes into account interaction. The non-interacting 
carriers are independent and fid — f^dip) + Since 
only the second derivatives of the HED are necessary, one 
can write fid{p) — pkT{\ 4- Inp) and fid{n) = nkT{l + 
Inn). 

The largest interaction term for the non-degenerate 
plasma describes correlation effect. It was calculated by 
Debye and HuckelQ in a form 



fc{n+p) = -(2eV3K'/')yVfcr(n+p)3/2. (19) 

This term is independent of the spectra of electrons and 
holes. In our approximation this is the only term that has 
non-zero d^f/{dndp) and contributes to the drag effect 
(See Eqs.dSlIl)). 

I also take into account the exchange interaction be- 
tween electrons in each ellipsoid, between heavy holes 
and between light holes. This interaction term has a 
higher power of T in the denominator of the HED than 
the correlation term. The thermodynamic potential den- 
sity ^{fi, T) for this interaction can be written in a form 
of the Wigner-Seitz integral (SeeQ) 



npinp2(fpi(fp2 
(Pl-P2)2(27r)6^ 



(20) 



where Up is the Fermi function that has the Boltzmann 
form in this case. To find the HED one should express 
chemical potentials through the density of carriers. 

For a conduction band consisting of g equivalent ellip- 
soids of rotation one gets 



fex {n) 



4:K,y^gm±kT' 



(21) 



where masses m± and TO|| , are perpendicular and parallel 
to the rotation axis of an ellipsiod, a = m\\/m±. At a > 1 



I{a) = 



1\pn arctan( yo— 1 



(22) 



For a parabolic valence band with light hole mi and 
heavy hole 



fex{p) 



-Ke nrp {m\ + mf) 



(23) 



My final result for corrections to a regular Einstein 
relation reads 



Q{n,T) = 1- 



e^h?nl{a) 



V2«;3/2(/cT)3/2 A^Km^g{kT)^ 



2K(fcT)2(m^/^ +m^/^)2' 



(24) 



The Dember field in terms of mobilities has a form 
1 ffkT e'J^ 



Ed 



{eU„ + eUp) VV " K3/2V2fcf^ 
TTe^n^iml + mf)Up e^h^I{a)Un 




KkT{rnfl^ + mf'^Y 2K^gm±kT 



Mobility of carriers in a plasma 

It is assumed above that mobilities of the carriers are 
controlled by the lattice scattering. But even in this case 
there are important corrections to these mobilities due to 
the interaction of electrons and holes. Each carrier is sur- 
rounded by a screening atmosphere of the opposite sign. 
This atmosphere is polarized by an applied electric field. 
The field of this polarization is opposite to the applied 
field so that the effective field acting on the carrier is less 
than applied field. It can be interpreted as a decrease 
of the mobility. The theory of this effect was created by 
Debye, Huckel, and Onsager (See Ref.[§]). The resulting 
changes of the mobilities are 



S{n,T) 



AC/„ A{/„ 



2^e3nV2 



3(1 



/0:5)At3/2(fcT)3/2 

(26) 

Coming back to the Einstein relation Eq. p7)) one should 
note that if the mobilities Un , Up are measured in the 
presence of plasma, the above corrections are irrelevant 
because the experimental values J7„,C/p contain them. 
However, if the mobilities are known from experiments 
without light excitation, as in the case of Ref. Q, the 
Einstein relation Eq. (fT7|) takes a form 



(27) 



e Up + Un 



where P{n,T) = Q{n,T) + S{n,T). One can see that 
this change increases the numerical factor in the second 
term of Q, originated from correlation, by 1.39. 



Discussion and Conclusion 

Now I discuss the low temperature anomaly of the 
CAD in Si. Function P{n, T) is shown in Fig.l in a proper 
temperature range. The comparison of the first approx- 
imation (correlation) with the second one (exchange) 
shows that the perturbation theory looks reasonable at 
T > 150K and at no = 1. 
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FIG. 1: (Color online)Function P{n,T) for Si at three differ- 
ent values of no defined as n = no x 2.3 x lO'^^cm"^; no = 2 
for lower curve (blue), no = 1 for two intermediate curves, 
and no = 0.5 for upper curve(green). The upper intermediate 
curve does not take into account exchange interaction. 

Da,cm^/sec 




FIG. 2: (Color online)Temperature dependence of Da ob- 
tained under assumption that Up ~ rp-2.2 _ 2 3 x 
lO^'^cm"^. The absolute values of the hole mobility is cho- 
sen such that maximum value of Da = 20cm^/sec. Maximum 
occurs at T ~ 230K, which is very close to the experimental 
result. 

To estimate CAD as a function of T one should know 
mobilities C/„ and Up. The experimental and theoreti- 
cal data of Ref.0 show that in silicon Up/Un ~ 0.25, at 
T « 200K. Then UpUn/{Up + f7„) « Up. The hole mo- 
bility in the pure silicon is due to the phonon scattering 
and it depends on temperature as T~^'^ in all tempera- 
ture range considered. The deviation from a usual law 
T~^-^ is due to the warping of the top of the valence 
band. The mobility of doped silicon with the hole den- 
sity 2 X has T-^-^ dependence^ at T > 180K 
that may be interpreted as a phonon scattering but with 
the warping smeared by the doping. A pure silicon is 
considered here, and r^2.2 mobility dependence is used. 

The expression Da = P{T)UpT with Up = RT-^ '^ 
is used to get T-dependence of Da that follows from 
the above theory. To make comparison with experi- 
mental result easier the factor R is chosen such that 
Da = 20cm^/sec in the maximum, similar to the experi- 
mental data of Ref . The theoretical result is shown in 
Fig. 2. 

Since the HED is calculated using perturbation theory. 



the discussion of the low temperature behavior might be 
doubtful, and the most important argument is position 
of the maximum of CAD. Clearly the way factor R is 
chosen has no effect on the maximum position. Theoret- 
ical position of maximum is 230K, which is close to the 
experimental position that has some uncertainty because 
of the large error bars. I think this similarity is a strong 
argument in favor of the proposed explanation. 

In the range T > 150K the the theoretical curve is sim- 
ilar to experimental points of Ref(0). At lower T theo- 
retical values become negative. This definitely means col- 
lapse of the perturbation theory, because negative CAD 
leads to an absolute instability of a neutral plasmaQ. It 
is intriguing that just near this temperature the character 
of experimental data changes: CAD becomes indepen- 
dent of both T and n. This might be a manifestation of 
a new phase. The speculations about this phase are out- 
side the scope of this paper. I would only mention that 
amount of excitons at these temperatures, as given by 
the Saha equation, is negligible but the Saha equation is 
not reliable for the case of non-ideal plasma. It is known 
also that the phase of exciton gas-liquid coexistence cor- 
responds to lower temperature at these densities (lo{. 

Finally, I proposed an explanation of the low temper- 
ature anomaly of the CAD based upon the Einstein re- 
lation for interacting carriers. The applicability of the 
theory is limited at low enough temperatures because 
the calculation of the HED is pcrturbational. Using fun- 
damental thermodynamics relations Eqs. p5ll7ll8p one 
could restore unknown thermodynamic functions of the 
non-ideal plasma at low temperatures by measuring mo- 
bilities, CAD, and the Dember field. 

I am grateful to M.I. Dyakonov for an important crit- 
ical comment and to Hui Zhao for a valuable discussion. 
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